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Abstract 



We present the first set of quenched QCD measurements using the 
recently parametrized fixed-point Dirac operator D FF . We also give a 
general and practical construction of covariant densities and conserved 
currents for chiral lattice actions. The measurements include (a) hadron 
spectroscopy, (b) corrections of small chiral deviations, (c) the renormal- 
ized quark condensate from finite-size scaling and, independently, spec- 
troscopy, (d) the topological susceptibility, (e) small eigenvalue distribu- 
tions and random matrix theory, and (f) local chirality of near-zero modes 
and instanton-dominance. 
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1 Introduction 



The only way to study a quantum field theory non-perturbatively from hrst 
principles is via lattice regularization. For a strongly-interacting theory like 
QCD, many essential features are non-perturbative. Over the course of more 
than two decades, lattice QCD has determined, with varying degrees of accu- 
racy, the hadron mass spectra, quark masses, the strong interaction coupling 
constant, low-energy constants, inter-particle potentials, the phase structure of 
the theory, decay constants, matrix elements and many other quantities M . To 
make lattice QCD into a precise science requires detailed study of the systematic 
error introduced by the lattice discretization. As observed recently in a study of 
the 2d 0(3) non-linear sigma model H, assumptions about how to extrapolate 
lattice results to the continuum may be wrong. Even for very accurate lattice 
measurements, large deviations as the lattice spacing a is varied could mean 
large uncertainty in continuum quantities, due to the extrapolation ansatz. To 
reduce such sensitivity, many groups now use e.g. C(a)-improved lattice actions 
to remove the leading order lattice artifacts. However, it may still be necessary 
to go to very fine resolution before the continuum extrapolation can be done 
confidently. 

A separate issue is that, for a long time, it was thought impossible to have 
chiral lattice fermions ||. This prevented a clean study of the chiral aspects 
of QCD, such as spontaneous symmetry breaking, and caused several technical 
headaches, such as fine-tuning the bare quark mass to make the pions light, 
exceptional configurations, mixing of operators with different chiral represen- 
tations and additional renormalization factors. The question was also raised 
whether a chiral theory, such as the Standard Model, could in principle be 
defined non-perturbatively. The domain wall fermions the related overlap 
construction j|] and the fixed-point (FP) action [|| are the known examples for 
lattice regularized fermions with chiral symmetry. The common feature behind 
these rather different constructions is the Ginsparg- Wilson (GW) relation [Q 
which is satisfied by all these constructions Q . Actually, one might construct a 
Dirac operator using directly the GW relation and the requirement of locality 
[[)) . The GW relation not only implies all the physical consequences of chiral 
symmetry |Io| , [Tlf , but also the existence of an exact symmetry transformation 
[12 . This has caused an explosion of interest in the last few years in theory 
[13 and in using different methods to explore the chiral behavior of quenched 
QCD (no sea quarks), which is now quickly maturing. With the standardly used 
algorithms, it is impossible to use chiral lattice actions for simulations of full 
QCD (with sea quarks), as these actions have a complicated structure |14fl . This 
is the major obstacle for the future use of chiral actions. For a recent promising 
development see, however, ref. Jl5[ . 

The fixed-point method is motivated by renormalization group properties of 
lattice field theories (l6|. Instead of trying to remove lattice artifacts order by 
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order in the lattice spacing a, a fixed-point (FP) action is designed to be less 
sensitive to the discretization even at large lattice spacings. In fact, for some 
quantities, a FP action is entirely blind to the discretization. An exact FP 
QCD action has exact chiral symmetry and offers a completely new approach to 
examine the chiral properties of the theory. Exact and approximate fixed-point 
actions have been constructed and tested for a number of models, including pure 
Yang-Mills theory, and in general have shown very good scaling behavior jl7|, [l8| . 
The construction and initial tests of an approximate fixed-point Dirac operator 
D FP are described in p9. Like other chiral actions, fixed-point actions are 
more costly to use in simulations than standard actions. However, FP actions 
are designed to have much reduced lattice artifacts and, if chiral actions are the 
correct long-term approach for lattice QCD, FP actions might be the optimal 
choice. 



In this paper, we present the results of a preliminary study of an approximate 
FP QCD action, to test if it is feasible to use such an action in simulations. A 
larger systematic study, which is running already as part of the program of the 
Bern-Graz-Regensburg Collaboration, is required to examine the lattice spacing 
dependence of measurements using this action (see ref. jl^]). As well as these 
initial measurements, we also give a practical method to construct and use 
conserved currents in simulations with chiral lattice actions. Many groups now 
use chiral actions, but non-conserved currents and/or non-covariant densities 
and so lose part of the advantages offered by a chiral formulation. 



This paper is a detailed description of the work summarized in 21, |22|. 
The paper is organised as follows. In Section 2, we examine how small chiral 
deviations of the action can be corrected. Hadron spectroscopy measurements 
with the parametrized FP Dirac operator are described in Section 3, giving the 
pion, rho and nucleon masses, the speed of light, the remnant additive quark 
mass renormalization and an indication of the presence of quenched chiral loga- 
rithms. In Section 4, we use the overlap-improved FP Dirac operator to measure 
the meson spectrum, the renormalized quark condensate via finite-size scaling, 
the quenched topological susceptibility, investigate random matrix theory and 
small eigenvalue distributions on small lattice volumes, and examine instanton- 
dominance of near-zero modes of the Dirac operator. The construction of covari- 
ant currents and densities and a discussion of Ward identities for chiral lattice 
actions are given in Section 5, followed by our conclusions. In the Appendix we 
collect some useful identities which are implied by the Ginsparg- Wilson relation. 
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2 Chiral behavior of D 



A Dirac operator D which satisfies exactly the Ginsparg- Wilson relation 

j 5 D + Dj 5 = D l5 2RD, (1) 

where R is any local function and trivial in Dirac space, has exact chiral symme- 
try even at non-zero lattice spacing a. Such a Dirac operator has infinitely many 



couplings 1 23 . Simulations are only feasible with a finite number of couplings, 
hence the chiral symmetry is only approximate. What needs to be examined is 
how good the approximation is. The exact fixed-point Dirac operator satisfies 
exactly the Ginsparg- Wilson relation. The parametrized fixed-point Dirac oper- 
ator D FP , which we use in simulations, is an approximate solution to the QCD 
fixed-point equations. The construction and initial tests of D FP are described 
in ||^, H^]. Here, we test the quality of the chiral symmetry of this operator. 

As will be discussed in Sec. 5, the function 2R can be absorbed into the def- 
inition of D, so for simplicity, we set 2R — 1 in most of the following equations. 
The Ginsparg- Wilson relation can be re-written as 

D + L>t = D ] D, (2) 

using the property 75-D75 — . If D satisfies this relation, its complex eigen- 
values lie on the circle of radius 1 and centre at (1,0) in the complex plane. 
In Fig. |l], we plot the complex eigenvalues of D FP for fixed-point gauge action 
S FP configurations. The lattice spacing can be determined via the Sommer 
parameter ro of the static quark anti-quark potential. To orient the reader, 
the bare coupling of the gauge action (3 — 6/g 2 — 3.0 and 3.2 correspond to 
lattice spacings a ~ 0.16 fm and a ~ 0.13 fm, respectively. The figure contains 
all the eigenvalues for 4 4 gauge configurations and those closest to the origin for 
8 4 configurations. It's clear that the breaking of the Ginsparg- Wilson relation 
is small over the entire eigenvalue spectrum even though the lattice spacing is 
quite coarse, and that the chiral symmetry remains very good even for larger 
volumes. For comparison, the eigenvalues of the Wilson Dirac operator D w at 
a similar lattice spacing are very far removed from the Ginsparg- Wilson circle 



Defining the operator A = 1 — D, the Ginsparg- Wilson relation is equivalent 
to the requirement that A be unitary i.e. A^ A = 1. We can measure how much 
this requirement is broken by finding the smallest and largest eigenvalues of 
A^A. In Fig. |[ we plot the ratio of the 50 smallest eigenvalues of A^A to 
the largest eigenvalue for a number of S FP gauge configurations with a lattice 
spacing 0.16 fm. The largest eigenvalue is approximately 1.5 and 38 for D FP 
and D w respectively. We see that the fixed-point operator is much closer than 
the Wilson operator to satisfying the unitarity condition. Also, as the volume 
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increases, it's more likely that using D produces a few very small eigenvalues, 
but the vast bulk of the A' A spectrum is close to 1. 

For some measurements, it's necessary that the remnant chiral symmetry 
breaking be very small. Given any suitable Dirac operator Dq as input, the 
overlap Dirac operator || 

D ov = 1 - A/VaTa, A = 1+s-D (3) 

is an exact solution of the Ginsparg- Wilson relation. The real parameter s can 
be used to optimize the convergence of approximations or the locality of the 
resulting overlap Dirac operator. The nice feature of this operator is that it's 
an explicit construction, unlike the fixed-point operator which is a solution to 
a set of equations. If the input operator Dq is already chiral, then A^A — 1 
and hence D ov — Dq. Most simulations using D ov have taken D w as the input 
operator, which has s evere ly broken chiral symmetry, as A' A is far from 1 and 
hence computing 1/vA^A is a very expensive numerical problem. A measure of 
this difficulty is the condition number of A* A, the ratio of its smallest to largest 
eigenvalue. Also, the overlap operator may inherit the large lattice artifacts of 
D w . Alternatively, taking D FP as the input operator, A^A is close to 1 and 
1 / V At A is easier to evaluate — the improvement in the condition number is 
clear in Fig. ||. The overlap construction should only bring small corrections to 
D FP . This way, the residual chiral symmetry breaking can be removed without 
destroying the very important fixed-point properties. 

As the chiral symmetry is already well approximated by D FP , we use a Leg- 
endre polynomial approximation of l/yAfA to construct the overlap-improved 
operator D FF . As A' A may occasionally have a few very small eigenvalues, 
we always project out the smallest 10 to 100 eigenvalues, which are treated 
exactly. One measure of the remnant chiral symmetry breaking is the vector 
norm Aq-w = \ \(D + D' — D'D)v\\, where v is a unit vector of random entries. 
If D solves the Ginsparg- Wilson relation, Aq-w vanishes. In Fig. ||, we plot 
Agw(A0 versus the order N of the polynomial used to approximate 1/v A^A 
for gauge configurations of volume 10 4 at different lattice spacings (0.16 fm, 
0.13 fm, 0.10 fm). The chiral symmetry breaking falls off exponentially as the 
polynomial order increases. The rate of the fall off varies depending on the 
range of the A^A eigenvalues which are not projected out. Projecting out the 
same number of A^A eigenvalues for the different lattice spacings, the rate of 
the fall off is the largest for the smallest lattice spacing. However, choosing the 
number of projected A^A eigenvalues such that ratio X m i n (A^ A) / X max (A^ A) for 
the remaining A^ A eigenvalues is approximately the same for the three differ- 
ent lattice spacings, the rate of the fall off is roughly equal for all the lattice 
spacings, which illustrates that the convergence rate is indeed governed by the 
ratio A m i n 

If the eigenvalues A n of D lie exactly on the Ginsparg- Wilson circle, then 
A„ = A„/(l — X n /2) are purely imaginary i.e. the eigenvalues are projected 
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from the circle onto the imaginary axis. In Fig. || , we plot |Re(A„)| versus the 
order of polynomial used to approximate for gauge configurations with 

lattice spacings 0.16, 0.13 and 0.10 fm. As the order increases, |Re(A„)| becomes 
exponentially small as the eigenvalues A„ lie closer and closer to the Ginsparg- 
Wilson circle, with the most rapid decrease at the finest lattice spacing. 

A lattice action must be local for its continuum predictions to be universal 
i.e. independent of the choice of lattice discretization. Locality on the lattice 
means that fields at large separation have an exponentially small coupling. The 
exponential decay of the action should be faster than that of the correlation 
functions. This issue was first addressed for the overlap operator D ov in fzEf , 
where a free parameter was varied to maximize the exponential fall-off. One 
measure of locality is 

/(r) = max{||D ov i;||, \\y - x\\ = r}, (4) 

y 

where v is vector with a point source at x and r is the square norm. In Fig. ^ we 
plot the expectation value of /(r)//(0) versus r for different overlap operators 
on 12 4 volumes at a lattice spacing of 0.16 fm. Compared with other tests of 
locality of D ov using D w as input (25J ||(| , the locality is significantly improved 
if the overlap operator is constructed using _D FP , with a faster exponential 
decay. The exponent v describing the decay at large separation r is v = 0.94 
for D^, whereas it is v = 1.60 for D FP . In order to optimize the locality 
of the parameter s in Eq. (|J) has to be tuned and we find s = 0.6 as 
optimal value at this lattice spacing. In contrast to there is no tuning of 
the parameter s needed for D FF to get the (approximately) optimal locality, 
i.e. we can use s = 0.0 for D FF . 

From these tests, we see that the chiral symmetry is well approximated by 
D FF and that the residual breaking can be removed using overlap-improvement. 
Only a relatively low order polynomial approximation is required to construct 
D FF , which is essential as D FP is more expensive to use than D w . The 
parametrized FP operator has 9 times more offsets and the computational cost 
per offset is also higher since D FP includes all elements of the Clifford algebra. 
If, in a matrix-vector multiplication, for D w the time needed for Dirac matrix 
multiplications is neglected, then D FP is a factor of « 4/offset more expensive 
leading to an estimated relative overhead of 36. Note, however, that this num- 
ber may vary quite a bit depending on the underlying computer architecture, 
as one of the main issues in present lattice simulations is rather fast memory 
access and fast communication than fast floating point units pTJ]. 
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3 Spectroscopy 



Measuring the mass spectrum is a standard benchmark for using a particular 
lattice action. Simulations with improved Wilson and staggered fcrmions in full 
QCD are now very advanced in reaching the physical theory. Mass spectroscopy 
in quenched QCD with chiral fcrmions is quickly maturing and is a very good 
test of the chiral behavior of the theory. In this preliminary study, we test the 
behavior of D FP , for example, how intact are the fixed-point properties, how 
large is the additive quark mass renormalization, how small an m w /m p ratio 
can be attained and do we observe quenched chiral logarithms. 



3.1 Simulation parameters 

The FP gauge action we use has previously been studied in detail To 
remind the reader, the bare coupling of the gauge action (3 = 6/g 2 — 3.0 and 

3.2 correspond to lattice spacings a ~ 0.16 fm and a ~ 0.13 fm respectively, 
as determined from the Sommer parameter ro of the static quark anti-quark 
potential [^8| . Lattice volumes of 6 3 x 16 and 9 3 x 24 are compared to investigate 
the volume dependence of zero mode effects. We use gaussian smeared sources 
and point sinks. Configurations were fixed to Landau gauge using the Los 
Alamos algorithm with stochastic overrelaxation p9|j . With the parametrized 
FP Dirac operator D FP , we simulate at input quark masses ma ranging from 
0.016 to 0.23. The smallest quark mass corresponds to m^/nip = 0.30(3). It 
has to be stressed that we introduce the mass in the most simple manner by 
defining^ D{m) = D FP + m or with the covariant scalar density in Eq. ( |S8| ) 
D{m) = (1 — m/2)D FP + m(2R)~ l , whereas the fixed-point Dirac operator for 
non-zero mass would in fact need a reparametrization by iteratively solving the 
RG transformations for every mass value as we did for the zero-mass Dirac 
operator D FP . Therefore our parametrization deviates more and more from 
the classical renormalized trajectory for larger masses. A multimass BiCGstab 
algorithm fl30| is used to invert the Dirac operator at all masses simultaneously. 
Errors are estimated with bootstrap resampling of the hadron correlators, which 
were symmetrized around t — T/2 to increase statistics. A correlated fit to 
the measured correlators is performed in the interval [i m i n , £ m ax] , where the 
minimum time i m ; n of the fit range is determined by two criteria: i) the effective 
mass starts to show a reasonable plateau, ii) the value of % 2 per degree of freedom 
in the fit (x 2 /df) starts to show a plateau and is of order 1. The maximum time 
t max is generally set to T/2 in the mass measurements, while for the finite 
momentum measurements it is reduced according to the length of the plateau. 

Fig. H shows the masses of the pseudoscalar and vector meson and the nu- 

2 The non-covariant form D FP +m for the massive Dirac operator was used for the results 
on the 6 3 X 16 lattice and the finite-momentum calculations on the 9 3 X 24 lattice 
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fermion action 


lattice size 


/3 


L s 


N con f 


ma 


parametrized FP 


6 3 x 16 


3.0 


1.0 fm 


96 


0.02,..., 0.20 


parametrized FP 


9 3 x 24 


3.0 


1.5 fm 


70 


0.016, ...,0.23 


overlap FP 


9 3 x 24 


3.0 


1.5 fm 


28 


0.012,..., 0.23 


overlap FP 


9 3 x 24 


3.2 


1.2 fm 


32 


0.012,..., 0.23 



Table 1: Simulation parameters for hadron spectroscopy. 



cleon measured with the parametrized FP Dirac operator on the 9 3 x 24 lattice. 
The inversion of the Dirac operator converged for all configurations within rea- 
sonable time, although for one configuration the number of inversions was about 
three times larger than the typical value of ~ 250. This indicates that the resid- 
ual chiral symmetry breaking leads to fluctuations on the order of our smallest 
quark mass for the real eigenvalues of the Dirac operator, and therefore it would 
not be wise to go to smaller masses than ma = 0.016 with the parametrized FP 
operator at this lattice spacing and volume. An extrapolation of the rho mass 
to the physical value gives a lattice spacing of a ~ 0.17 fm, which is slightly 
higher than that obtained from the Sommer parameter. The sign and size of 
this deviation is consistent with earlier findings J3l| . Fig. [?] is an Edinburgh plot 
from this data. 



3.2 The pion mass in the chiral limit 

A lattice calculation of the pion mass in the limit of small quark mass with the 
FP action is interesting for several reasons: first, it serves as a quantitative check 
of how chiral our action actually is. For an exactly chiral fermion action like 
the overlap or the exact fixed-point action, there is no additive renormalization 
of the quark mass. But does the approximation we make by parametrizing 
the FP action introduce a residual mass, and can we measure it? Second, at 
the lattice volumes we are working with, one expects to see topological finite- 
volume effects proportional to l/(m^W) and l/(m q W) |32) which complicate 
a reliable pion mass measurement at small quark masses. It is therefore crucial 
to check whether these effects are under control. Third, having an action with 
good chiral properties at hand, one can check for the appearance of logarithmic 
terms in the chiral extrapolation of the squared pion mass. 

We use three different operators to extract the pion mass: the pseudoscalar 
P = ^75 the fourth component of the axial vector = ^7475 \t and the 
difference of pseudoscalar and scalar S = It'll/. In the quenched theory, the 
pseudoscalar correlator (P(x)P(0)) is contaminated by the topological finite 
size effects from configurations with non-trivial topology, as mentioned above. 
The same special finite size contributions appear in the scalar correlator, so by 
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ma 


m 7r (P) 


X z /df 


mn (A) 


X 2 /df 


m w (P-S) 


X 2 /df 


0.02 


0.394(41) 


0.77 


0.345(62) 


1.55 


0.273(33) 


1.19 


0.03 


0.398(25) 


1.01 


0.381(30) 


0.62 


0.341(32) 


0.72 


0.04 


0.416(19) 


0.92 


0.403(21) 


0.64 


0.381(28) 


0.47 


0.05 


0.439(16) 


0.87 


0.428(19) 


0.80 


0.418(22) 


0.33 


0.06 


0.463(14) 


0.90 


0.451(17) 


0.98 


0.451(21) 


0.25 


0.08 


0.502(15) 


0.91 


0.498(14) 


1.14 


0.512(18) 


0.19 


0.10 


0.546(13) 


0.74 


0.541(12) 


1.10 


0.566(16) 


0.22 


0.13 


0.611(12) 


0.53 


0.603(10) 


0.93 


0.638(16) 


0.37 


0.16 


0.672(11) 


0.38 


0.663(10) 


0.78 


0.703(15) 


0.64 


0.20 


0.748(10) 


0.26 


0.737(08) 


0.64 


0.783(13) 


1.11 



Table 2: Pion masses for parametrized FP fermions on 6 3 x 16 lattice from three 
different pion operators. The quark mass is introduced by D FP (m) = D FP + re- 
building the difference of the pseudoscalar and scalar correlator (P(x)P(O)) — 
(S(x)S(0)}, this artificial effect should cancel, although if the action is not 
exactly chiral, the cancellation is also not exact. For larger quark masses, the 
P-S correlator is contaminated by effects from the scalar part and is expected 
to deviate from the pion. For the axial vector correlator (A4(z)Ai(0)), the 
divergent contribution is expected to be partially suppressed J33, M . 



ma 


77V (P) 


X 2 /df 


m n (A) 


X 2 /df 


777, (P-S) 


X 2 /df 


0.016 


0.267(20) 


1.25 


0.223(15) 


0.53 


0.209(21) 


0.64 


0.02 


0.274(15) 


1.18 


0.256(14) 


0.72 


0.243(16) 


0.66 


0.025 


0.292(11) 


0.92 


0.281(10) 


0.68 


0.272(13) 


0.72 


0.03 


0.314(8) 


0.93 


0.306(9) 


0.88 


0.298(11) 


0.80 


0.04 


0.351(7) 


0.98 


0.346(8) 


1.12 


0.342(9) 


0.89 


0.05 


0.384(6) 


1.06 


0.381(7) 


1.23 


0.379(8) 


0.97 


0.06 


0.415(6) 


1.13 


0.413(7) 


1.25 


0.413(7) 


1.06 


0.08 


0.473(6) 


1.27 


0.473(6) 


1.20 


0.476(6) 


1.24 


0.10 


0.526(5) 


1.38 


0.526(6) 


1.16 


0.532(5) 


1.41 


0.13 


0.601(4) 


1.52 


0.600(5) 


1.16 


0.609(5) 


1.69 


0.17 


0.692(4) 


1.71 


0.691(5) 


1.30 


0.703(4) 


2.00 


0.23 


0.819(3) 


2.00 


0.817(4) 


1.76 


0.832(4) 


2.16 



Table 3: Pion masses for parametrized FP fermions on 9 3 x 24 lattice from 
three different pion operators. 

Figs. ^ and ^| show a comparison of the squared pion mass versus the input 
quark mass at the two lattice sizes 6 3 x 16 and 9 3 x 24 using the parametrized 
FP Dirac operator. On the smaller volume, the three different pion correlators 
give very different values at small quark masses. The P correlator lies highest, 
while the P-S correlator gives considerably smaller pion masses, as the inset 
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in Fig. g shows. The A correlator lies in between the two. This behaviour is 
in qualitative agreement with the results from domain- wall fermions |32[ . For 
larger quark masses, the P-S correlator deviates from the other two, as it is then 
difficult to disentangle the contribution of the closely lying scalar state in the 
mass determination, which can also be seen from the increasing value of x 2 f° r 
the mass fit in Table [| As expected, the discrepancy between the three corre- 
lators in the chiral limit is reduced at our larger lattice volume. The ordering 
among the correlators is still the same, but the pion mass difference is smaller. 
A recent study with Wilson overlap fermions at a similar spatial lattice size of 
« 1.4 fm also shows no significant difference |j6|, but their smallest quark mass 
is much larger than ours. We can conclude from our findings that at the rather 
small statistics that we have collected, the topological finite-volume effects be- 
come non-negligible at the smallest few quark masses, where the pseudoscalar 
correlator is clearly contaminated and cannot be used to get reliable results. 
We therefore decide to work in the following with the P-S correlator at small 
and with the P correlator at large quark mass, changing correlators at an inter- 
mediate quark mass where both agree, which is at ma = 0.06 for the 9 3 x 24 
lattice. 



ma 


m n /m p 


m p 


X 2 /df 


m N 


X 2 /df 


0.016 


0.304(33) 


0.673(21) 


0.87 


0.919(117) 


0.66 


0.02 


0.357(26) 


0.679(21) 


0.48 


0.928(102) 


0.83 


0.025 


0.396(20) 


0.686(14) 


1.06 


0.908(74) 


0.95 


0.03 


0.429(17) 


0.693(12) 


1.42 


0.927(55) 


1.03 


0.04 


0.482(15) 


0.708(12) 


1.79 


0.975(43) 


1.43 


0.05 


0.523(14) 


0.723(11) 


2.02 


1.013(32) 


1.82 


0.06 


0.562(11) 


0.735(10) 


2.20 


1.043(27) 


2.38 


0.08 


0.619(10) 


0.763(9) 


2.45 


1.098(19) 


2.97 


0.10 


0.666(9) 


0.788(9) 


2.56 


1.152(15) 


3.22 


0.13 


0.726(8) 


0.828(7) 


2.51 


1.232(11) 


3.46 


0.17 


0.786(7) 


0.886(7) 


2.28 


1.335(11) 


3.71 


0.23 


0.852(6) 


0.979(6) 


2.16 


1.489(10) 


4.17 



Table 4: m^/mp, rho and nucleon masses for parametrized FP fermions on 
9 3 x 24 lattice. 

On the larger lattice, we also measure the unnormalized AWI quark mass 

qi ) - EAP^mo)) ' (5) 

and take the average of the ratio of correlators over the range 5 < t < 9. We use 
an ultralocal (non-conserved) axial current and neglect the renormalization 
factors Za and Zp which would show up on the right hand side of Eq. (|^), as 
they are not relevant for the following analysis. We fit the data with a linear 
function to determine whether the remaining chiral symmetry breaking of the 
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action introduces a residual mass. The smallest two masses were left out of the 
fit, which is shown as a dashed line in Fig. |9|. From the intersection of the linear 
fit with the horizontal axis we read off a residual mass m^ l a = —0.003(4), 
which is consistent with zero within errors. A quadratic fit to the squared pion 
mass from the P-S correlator for am < 0.05 and the P correlator for am > 0.05 
gives m res a = 0.004(4) with a value of \ 2 /df = 0.28, while a fit to the form 
predicted by quenched chiral perturbation theory (Q^PT) p9| 



to ps = 2A(m q + m res ) 



1 



/ 2A(m q + m rcs ) \ 
V A 2 + J 



4B(m q + m rcs ) 2 , 
(6) 



gives m les a = —0.004(5) and 0.23(7) < S < 0.30(18) when the scale A x is varied 
in the range between 0.8 GeV and 1.2 GeV, with typical values of x 2 /df = 0.05. 
Obviously the large errors do not allow to single out a preferred form for the 
chiral fit, but for the QxPT form, the residual quark mass agrees with the one 
from the axial Ward identity, while the agreement is worse in the case of the 
quadratic fit. This gives a hint that we see a signal of the chiral logarithm in our 
pion mass measurements. (Note that the smallness of the x 2 /df values is due 
to the fact that the results for different quark masses are strongly correlated, 
hence an improvement from 0.28 down to 0.05 is meaningful.) 



3.3 Energy-momentum dispersion relation 

As the FP action is classically perfect, its dispersion relation is expected to 
show small scaling violations. In Fig. [l(] the squared speed of light c 2 = 
(E(p) 2 — m 2 )/p 2 measured on the 9 3 x 24 lattice is shown for the smallest 
momentum p = 27rn/9 with \n\ = 1 (|p| ~ 0.8 GeV) at different quark masses. 
While c 2 for the pion is consistent with 1 within errors, it is not for the rho 
meson. However, the (purely statistical) error bars do not include systematic 
uncertainties from choosing the fit range, which was especially difficult for the 
rho in this low statistics, small volume study. As a qualitative result, it is clear 



that compared to the Wilson or Sheikholeslami-Wohlert clover action 37 1, the 
dispersion relation is significantly improved. 



4 First results for the overlap- improved FP ac- 
tion 

Using the parametrized FP Dirac operator D FP as a starting point for an overlap 
expansion, it is possible to remove any effects from residual chiral symmetry 
breaking which is due to the imperfection of the parametrization. Note that 
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ma 


n / Li 


m 7r (P) 


run (A) 


mJP-S) 


m p 


0.012 


0.273(38) 


0.235(12) 


0.232(15) 


0.202(23) 


0.741(56) 


0.016 


0.313(32) 


0.257(19) 


0.258(13) 


0.234(20) 


0.747(43) 


0.02 


0.348(29) 


0.279(10) 


0.280(12) 


0.261(18) 


0.752(36) 


0.03 


0.418(24) 


0.326(8) 


0.328(11) 


0.319(14) 


0.763(27) 


0.04 


0.472(22) 


0.370(7) 


0.370(9) 


0.365(11) 


0.773(26) 


0.06 


0.556(18) 


0.444(6) 


0.443(7) 


0.443(8) 


0.797(22) 


0.09 


0.644(15) 


0.536(6) 


0.534(6) 


0.541(5) 


0.832(17) 


0.13 


0.723(12) 


0.640(5) 


0.639(6) 


0.651(5) 


0.886(13) 


0.17 


0.776(10) 


0.734(4) 


0.733(5) 


0.748(5) 


0.946(10) 


0.23 


0.829(7) 


0.863(4) 


0.862(5) 


0.880(4) 


1.042(8) 



Table 5: to w /to p , pion and rho masses for overlap- improved FP fermions on 
9 3 x 24 lattice at (3 = 3.0 (lattice spacing 0.16 fm). 

the true FP Dirac operator would reproduce itself by the overlap procedure. 
Since our D FP is close to the true FP Dirac operator we expect that the overlap 
procedure will not drive it too much away from the FP, i.e. by improving chiral 
properties we don't spoil the good scaling properties. We define the massive 
overlap Dirac operator (for a discussion of the 2R ^ 1 case, see Sect. 5) 

D ov (m) = (l-^)D ov (0) + ^, (7) 

where the massless overlap operator is 

1 / A \ 1 
D ov = -= 1 - -== 8 
y/2R \ VaHJ VM 

and in the kernel, 

A = 1 - V2RD FP V2R, (9) 

the parametrized FP Dirac operator enters. R is an explicit parametrization of 
the fixed-point R in the Ginsparg- Wilson relation in Eq. (jl]). 

4.1 Meson spectrum 



For spectroscopy we approximate the inverse square root 1 / v A^A with a polyno- 
mial of order 3, projecting out the smallest 20 eigenvalues of A^A for faster con- 
vergence. The improved chiral behaviour allows the quark mass to be decreased 
even further. The smallest input quark mass we considered is ma — 0.012 where 
the pion to rho mass ratio is m^/mp — 0.27(4) and 0.24(7) at lattice spacings 
0.16 and 0.13 fm respectively (Tables § and g). In Figs. |ll| and |l|, we show 



12 



ma 


n / IJ 


m 7r (P) 


m v (A) 


mJP-S) 


m p 


0.012 


0.236(69) 


0.216(16) 


0.189(27) 


0.145(37) 


0.614(83) 


0.016 


0.297(61) 


0.232(14) 


0.214(22) 


0.176(30) 


0.592(66) 


0.02 


0.346(52) 


0.250(14) 


0.236(21) 


0.203(35) 


0.588(53) 


0.03 


0.446(41) 


0.296(11) 


0.287(19) 


0.264(18) 


0.592(36) 


0.04 


0.518(35) 


0.336(9) 


0.329(19) 


0.314(15) 


0.606(29) 


0.06 


0.624(28) 


0.406(9) 


0.405(16) 


0.398(11) 


0.639(23) 


0.09 


0.714(21) 


0.493(8) 


0.497(12) 


0.497(9) 


0.691(17) 


0.13 


0.782(17) 


0.593(7) 


0.596(8) 


0.605(7) 


0.758(14) 


0.17 


0.829(13) 


0.684(6) 


0.686(6) 


0.699(7) 


0.826(11) 


0.23 


0.873(10) 


0.811(6) 


0.811(5) 


0.828(6) 


0.929(9) 



Table 6: to w /to p , pion and rho masses for overlap- improved FP fermions on 
9 3 x 24 lattice at (3 = 3.2 (lattice spacing 0.13 fm). 



the pion mass squared measured using the overlap-improved FP Dirac opera- 
tor. For an exactly chiral action, the topological finite- volume effects cancel 
exactly for the P-S correlator. A quadratic fit to the squared pion mass is con- 
sistent with zero at m = within the large error at this small number of gauge 
configurations. 



4.2 Chiral condensate 



It is the common expectation that for QCD with a number Nf > 2 of mass- 
less quark flavours, the chiral symmetry is spontaneously broken by a non-zero 
expectation value for the chiral condensate (^$). Chiral perturbation theory 
(xPT), which is based on this assumption, is an excellent description of many 
low-energy QCD phenomena |Q. However, it is only possible via lattice QCD 
to test from first principles if the symmetry is spontaneously broken. 

The leading order effective theory of xPT contains the low energy constants 
/tt and E. In full QCD the Gell-Mann-Oakes-Renner relation 

flml = 4mS, (10) 

becomes exact in the m — > chiral limit and the chiral condensate (\I f \f r ) = 
(in) = (dd) = . . . , defined at zero quark mass, is equal to The low energy 
constant £ depends on the number of massless flavours Nf. 

In quenched QCD (Nf = 0), the relation (f*) = -E and Eq. (juj receive 
corrections even in the chiral limit |39, [h|. Actually, ('5^') and m^Jm in the 



chiral limit are not defined due to diverging quenched chiral logarithms. This 
expectation seems to be confirmed in numerical studies of the Banks-Casher 
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relation Q. 



On the other hand, it is possible to study and determine the low energy 
constant E(iV/ = 0) in the quenched theory. Under the assumption that S(iV/) 
is a smooth function of Nt and E(iV/ = 0) is close to E(iV/ = 3) (which is 
the standard assumption when quenched results are used to estimate full QCD 
quantities) we get an estimate for the chiral condensate — (WW) (Afy = 3) = 
E(N f = 3) « H(N f = 0). 

One possibility to determine E(./V/ = 0) is to study the chiral condensate 
in a fixed topological sector with charge Q in a finite volume V at finite quark 



mass m |40 . The volume and the quark mass are chosen so that the finite 
size effects are dominated by the pions with zero momentum. Using xPT, or 
random matrix theory, the fermion condensate at finite volume and quark mass 
has been calculated in the continuum, both for full and quenched QCD. The 
quenched QCD condensate (WW) mi y.Q is given by 

-<**W,Q = mUS 2 [/| Q |(z)^| Q |(z) + I lQl+1 (z)K lQhl (z)} + M , (n) 

where Iq and Kq are modified Bessel functions, z — mEU and the Nf = 0, 
V = oo, m = low energy constant E is the quantity we wish to measure. 
By measuring (WW) m ,v.Q m different topological sectors at different masses and 
volumes, the continuum prediction of the m and V dependence can be used to 
extract E. 

On the lattice we calculate the bare subtracted condensate by measuring the 
trace 



(D{m)2R)- 1 - -D(0)D~ x (m) 



' ' 'Tr' 



VI- f 



P(m)2i?)- 1 - i] ) . (l->) 



where D and i? are related via the Ginsparg- Wilson relation and where we make 
use of Eq. (Q) to simplify the expression on the first line. Notice that in order 
to keep the conventional notation we use {^^)m h vQ in Eq. ( |l2| ) even though 
we actually measure the expectation value of the scalar density S° as defined 
in Eq. (|53]), which differs from (^^)m b vQ by a tree level factor that will finally 
drop out in the renormalized value for E. We measure the trace stochastically 
using random Z(2) vectors In order to measure the condensate at very 

small quark mass, the remnant explicit chiral symmetry breaking must be very 
small, so we use the overlap- improved operator , where D FP is taken as 
input. Due to the \Q\ chiral zero modes of the Dirac operator, the quenched 
condensate contains a term |Q|/(mV) which diverges as the mass becomes small. 
In the volumes we study, the \Q\ zero modes always have the same chirality and 
their contribution to the condensate is removed by measuring the trace Tr' in 
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the chiral sector opposite to the zero modes [^3|, i.e. if the \Q\ modes have 
chirality +, the Z(2) vectors used to measure the trace are chosen to have 
chirality — . To determine ( , J n J / )™ 1 y q, the stochastic trace is doubled to include 
both chiral sectors. 

We first generated ensembles of gauge configurations at different volumes. 
We determined the topological charge Q of the configurations from the chirality 
of the lowest eigenmodes of (D^yD^y, which we find with an Arnoldi solver 
|l5). To scan for the topology, a Legendre polynomial of order 2 is used to 
approximate 1/VaTA, projecting out the 10 smallest A^A eigenvalues. This is 
sufficiently accurate to determine the chirality of the eigenmodes to better than 
1% accuracy. 



L 


order 


101 


^conf 


8 


7 


1 


154 






2 


41 


10 


10 


1 


53 






2 


43 



Table 7: Polynomial order of 1/VaTA approximation and statistics for the 
condensate measurement. 

We have measured the condensate in volumes 8 4 and 10 4 at lattice spac- 
ing 0.13 fm (corresponding to the bare coupling (3 — 3.2). We use 10 ran- 
dom Z(2) vectors to measure the trace for each configuration and a BiCGstab 
multi-mass solver to invert at all masses simultaneously. In the overlap 
operator, we approximate 1/v A^A with Legendre polynomials of order 7 and 
10 for volumes 8 4 and 10 4 , respectively. This gives sufficiently precise chiral 
symmetry — increasing the polynomial order further, the relative change in 
((**)™ b VQ a 3 )/(ma) is < 0(1O~ 4 ). We project out the 10 smallest A 1 A eigen- 
values, which are treated exactly. Our statistics are given in Table 0. In Fig. [l3], 
we plot (— (^ty^y Qa 3 ) / (ma) as a function of ma for the different volumes and 
topological sectors. 

The bare quark condensate at finite quark mass contains a ~ m/a 2 cut-off 
effect. As the quark mass m — * 0, 

m 2|Q| a 2 ' { ' 

where ci is an unknown coefficient which has to be fitted. As the coefficient c\ 
comes from ultraviolet fluctuations, it is natural to assume that it is independent 
of the topological charge Q. The contribution of the \Q\ zero modes has been 
removed, and there is no artifact l/(ma 3 ) in Eq. ( |l3| ) due to the fact that 
the condensate is defined as the expectation value of a scalar operator which 
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transforms covariantly under chiral transformations and has no mixing with 
the unit operator. (For further discussions, see Sec. 5.) From Fig. O, we 



see that (-(WW) 



sub 

m,V,Q 



'')/{r 



reaches a plateau at very small quark mass. 



In Fig. |I3, we plot the value of (-(WW)^" b vo a 3 )/(ma) at ma = 10~ 4 versus 



Vy(2|Q|), which we fit to the form of Eq. (|13|). From the slope, we extract 
the bare low energy constant as a 3 S = 4.42(36) x 10~ 3 . We convert this using 
the Sommer parameter (fa/a = 3.943(60) has previously been measured at this 
lattice spacing), giving TqE = 0.271(22)(12), where the first error is statistical 
and the second the uncertainty in the scale. 



In order to turn this bare result into the renormalized low energy constant 
we need the scalar renormalization factor Zg. In the present test study we ob- 
tained Zg combining the continuum extrapolated renormalization group invari- 
ant (RGI) quark mass of the ALPHA collaboration 0] with our spectroscopy 
data in Sect. 4.1. This method has been suggested recently by Hernandez et 
al. El. 



Chiral symmetry connects the scalar and pseudoscalar renormalization fac- 
tors: Zp = Zg = 1/Z m , where Z m is the multiplicative mass renormalization. 
This way the problem is reduced to finding Z m . The technique requires the 
measurement of the pion mass at a number of quark masses. The renormal- 
ization factor Z m connecting our quark mass m to the RGI mass M defined at 
some reference pion mass m w is given by 

Um 



(14) 

where Um = tqM . The renormalization factor Z m should be independent of 
the reference point, so finding this ratio at a number of pseudoscalar masses 
indicates the systematic error. We have performed mass measurements at two 
lattice spacings corresponding to j3 = 3.0 and 3.2. We use the same reference 
points and values of Um &s in [P^l to determine Z m . The results are summarized 
below in Table pL 



& 


*£rcf 


U M 


TqTTI 




3.0 


1.5736 

3.0 

5.0 


0.181(6) 
0.349(9) 
0.580(12) 


0.156(4) 
0.296(8) 
0.490(14) 


1.16(5) 
1.18(5) 
1.18(4) 


3.2 


1.5736 

3.0 

5.0 


0.181(6) 
0.349(9) 
0.580(12) 


0.141(27) 
0.271(24) 
0.446(22) 


1.28(25) 
1.29(12) 
1.30(7) 



Table 8: The renormalization factor Z m determined at different lattice spacings 
and reference pseudoscalar masses. 

At each lattice spacing, the values for Z m at the different reference points 
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are in very good agreement with one another. We take the average of the values 
and the error at a; re f = 3.0 as our determination of Z m . The renormalization 
group invariant low energy constant is E = ZsE. To convert this result into the 
MS scheme, we use the fact that 



m m (2GcV) 



E M s(2GeV) M 



0.72076, (15) 



where the ratio of masses has been calculated perturbatively to four loops in 
the MS scheme @. 

Taking the renormalization factor Z$ = 1/Z m = 0.78(7) at (3 — 3.2, where 
the bare E is obtained, we get for the renormalization group invariant low energy 
constant t-qE = 0.210(17)(9)(20), with statistical, scale and renormalization 
errors respectively. Taking r$ — 0.5 fm and combining the errors in quadrature, 
this corresponds to E = (235 ± 11 MeV) 3 . We can also convert this result to the 
MS scheme giving E^g(2 GeV) = (262 ± 12 MeV) 3 . As Fig. [TJ illustrates, our 
measurement of the renormalized chiral condensate is in good agreement with 
other recent determinations (3^, 44, 5l|. 



4.3 Alternative determinations of the chiral condensate 

A direct determination of E discussed above is only possible with a chirally sym- 
metric action. An alternative, more phenomenological way is to use the GMOR 
relation, Eq. (M) using the pion mass measurements and the experimental value 
°f fn fjsl, [361 . Although the ratio m^/m q is not defined in the quenched 
theory in the chiral limit m q — > 0, due to quenched chiral logarithms, we see 
that the pion mass squared is consistent with linear behavior for intermediate 
quark mass. Assuming that the chiral logarithms have very little effect in this 
mass range, we identify the slope Bm as 4E//^. Table || gives the value of the 
slope Bm and the bare condensate a 3 E calculated through the GMOR relation, 
E = B M fl/A, using the physical value f v = 131 MeV. The value at /3 = 3.2, 
a 3 E = 5.8(2) x 10 -3 has to be compared with the direct measurement discussed 
above, which gave 4.4(4) x 10~ 3 . 



/3 


a[fm] 


B M 


a 3 E x 10 3 


3.0 


0.16 


3.22(12) 


8.8(3) 


3.2 


0.13 


2.90(12) 


5.8(2) 



Table 9: The value of the bare chiral condensate determined from pion mass 
measurements. The slope Bm of rn^ vs. m q , and the value of the bare chiral 
condensate determined from GMOR relation, with = 131 MeV. 

We have to mention a difficulty with the technique used in the direct deter- 
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urination of the condensate discussed in the previous subsection. Finding the 
finite- volume behavior of the trace in Eq. (|l2|), the major problem is from rare, 
but very large contributions to the average. These come from configurations 
where the Dirac operator has very small, but non-zero, eigenvalues. Due to 
their presence, we found it difficult to control the statistical error. The distri- 
bution of the low-lying eigenvalues of the Dirac operator offers an alternative 
method to determine X | p2[ . This method does not suffer from the problem 
mentioned above. The appearance of very small eigenvalues does not bring any 
large contributions, which might make this method competitive, or even better 
than the finite-volume technique. 

To compare techniques of measuring the low energy constant S, we per- 
formed a test study of the distribution of low-lying eigenvalues of the Dirac 
operator. We generated ensembles of 4 4 lattice volume, with 2000 gauge config- 
urations at (3 = 2.4 (lattice spacing 0.30 fm) and 1000 configurations at /3 = 2.7 
(lattice spacing 0.22 fm). Admittedly, this lattice volume is too small and the 
resolutions are too coarse to allow a serious study of the low lying eigenvalues. 
Using the Ritz functional method, we measured the smallest non-zero eigenvalue 
of and determined the smallest eigenvalue distribution for each topological 
sector. For the ove rlap- improved operator, we used an order 5 polynomial ap- 
proximation of l/V A^A, with the smallest 20 and 8 A^A eigenvalues projected 
out at (3 = 2.4 and 2.7 respectively. Even on these coarse lattices, this low-order 
approximation has very little chiral symmetry breaking. Using random matrix 
theory, the distribution of the smallest eigenvalue A in topological sector Q in 
the quenched theory is 

P{z) = V 22 / 4 det[/ 2+ ^(z)] z,j = l,...,|Q| 

P(z) = ^e^ 2 / 4 Q = 0, (16) 

where z = AEU is the rescaled eigenvalue and Ik are the modified Bessel func- 
tions. In Fig. |l6|, we plot the measured distributions for topological sectors 
|Q| = 0,1,2 at (3 = 2.4, as well as the fits to the predicted form. We see 
the fits describe the data quite well, allowing us to estimate the bare quantity 
a 3 S = 0.0687(27), or alternatively r^E = 0.332(16). In Fig. [I?], however, we 
see that the smallest eigenvalue distributions at (3 — 2.7 for different topological 
sectors lie on top of one another and are clearly not described by Eq. ([!(]). The 
extreme values of (3 and/or lattice size (L — 4) might be the reason for this 
result. Further studies are needed to clarify this point. 



4.4 The quenched topological susceptibility 

As discussed above, in the course of measuring the quark condensate we have de- 
termined the topological charge of the gauge configurations from the chirality of 
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the lowest eigenmodes of (D^V ■ As a byproduct, we obtain the quenched 
topological susceptibility xt = (Q 2 ) /V from the distribution of the topological 
charge. From an ensemble of 200 10 4 configurations at our smallest lattice spac- 
ing 0.13 fm, we find r\xt — 0.0612(75), corresponding to xt = (196 ± 6MeV) 4 . 
This result, which still might have sizeable cut-off effects, is consistent with 
earlier determinations, as shown in fig. CM 



4.5 Local chirality of near-zero modes 



Exact zero modes of the Dirac operator tell us, via the index theorem, about 
the topological charge Q of the background gauge configuration. However, the 
exact zero modes alone cannot break chiral symmetry spontaneously. According 
to the Banks-Casher relation p3| ( 1 I r \f r ) = —7rp(0) ^ 0, the Dirac operator must 
build up a finite density of near-zero modes, which does not vanish as V — > oo. 

One mechanism which explains the formation of near-zero modes involves 
instantons. Consider a gauge configuration containing one instanton and one 
anti-instanton. If the instanton and anti-instanton are separated by a large 
distance, the Dirac operator has a pair of complex eigenvalues lying close to 
0. The farther the instanton and anti-instanton are separated from each other, 
the closer the complex eigenvalue pair moves to the origin. If the instanton and 
anti-instanton are brought closer together, the complex eigenvalue pair moves 
away from the origin and disappears into the bulk of the eigenvalue spectrum. 
If the gauge configurations contain many instantons and anti-instantons, this 
could produce a non-zero density of near-zero modes, giving p(0) ^ in the 
infinite volume limit. 



The question has recently been raised if it is possible to show that the near- 
zero modes are dominated by instantons. From instanton physics, it is expected 
that the modes are highly localized where the instantons and anti-instantons sit. 
If this is so, then in these regions the modes should be close to chiral i.e. mostly 
either left- or right-handed, depending on whether it is sitting on an instanton 
or anti-instanton. In [ |54| , the authors defined a measure of local chirality at 
lattice site x by 



tan 



-(1 + X(x)) 



*R f *R(a;) 



(17) 



where ^l/r^) are the standard L/R projections of the corresponding wave 
function. An exact zero mode is purely either left- or right-handed, giving 
X(x) = ±1 at all lattice sites x. If a near-zero mode is localized around 
instanton-anti-instanton lumps, then X(x) should be close to ±1 for the sites x 
where the probability density ^^{x) is largest. 
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In the original paper [ p4| . many near-zero modes of the Wilson operator 
D w were analysed for many gauge configurations and the finding was that in 
the regions where the modes are localized, the distribution for X(x) is peaked 
around and the modes do not display local chirality. This led to the conclusion 
that the near-zero modes are not dominated by instantons. Since then, several 
other groups have found the opposite conclusion |55| , using a Dirac operator with 
much better chiral symmetry than Z) w (or even just an alternative definition of 
a complete basis for the non-normal operator Z? w ). They find the distribution 
of X(x) is double-peaked with maxima at large positive and negative values of 
X, indicating that the modes are locally chiral. On the other hand, in several 
more detailed comparisons the semiclassical expectations were not confirmed by 
the numerical dataQ. Presumably one must not take this picture too seriously 
in a situation where these objects do not form a dilute gas and live in a strongly 
fluctuating background. 



We have analyzed the 10 smallest near-zero modes of the overlap-improved 
for 60 10 4 gauge configurations at lattice spacing 0.13 fm. We use a Legen- 
dre polynomial of order 2 to approximate 1/v A* A in , with the 10 smallest 
A^A being projected out and treated exactly. The eigenvalues and eigenvectors 
are found using the Arnoldi solver. In Fig. [l9| we plot the distribution P(X) of 
the measure of local chirality X at the lattice sites where the density ^^(x) of 
a mode is largest. The three distributions correspond to taking, for each mode, 
1%, 5% and 10% of all lattice sites which have the largest density ^^(x). We 
do not include exact zero modes, for which X(x) — ±1 at all lattice sites. We 
see a very clear double-peaked distribution, whose maxima are farther from zero 
if we only include the sites where the modes are most localized. The maxima are 
not at X = ±1 as the modes are not exactly chiral. We find the same conclusion 
as [[35] — where the near-zero modes are most localized, they are also very chi- 
ral. Recent work 1 57 has shown that at the places where the near-zero modes of 
the Dirac operator are concentrated, the gauge field appears to contain lumps 
of FF. This behavior is consistent with the picture of instanton-dominance of 
the near-zero modes, however it is not a conclusive evidence that instantons are 
the driving mechanism for chiral symmetry breaking. 



5 Covariant densities and conserved currents 

The GW relation^ 

y 5 D + D l5 = D l5 D . (18) 



3 The case 2R ^ 1 in Eq. (hi) will be considered at the end of this section 
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implies an exact SU(iV/) x SU(Nf) global symmetry on the lattice |Q. The 
vector transformation reads 

S^if, = ieT a i) , 8$$ = -i^T a t , (19) 

while the axial transformation has the form 

S a A ^ = ieT a %^ , 5%$ = UfosT** , (20) 

where T a , a = 1, . . . , AT? — 1 are SU(iV/) generators 

[T a , T b ] = if abc T c , tr (T a T b ) = U ab , (21) 

and 

75 =75(l-£>). (22) 

The action 

A = ^D(U)ip (23) 

is invariant under these transformations if D(U) satisfies Eq. ([l8|). The fact that 
A is a scalar under the transformations in Eqs. (jlj]) and ( p0| ) implies also that 
it is 0{a) improved since the mixing of the action density with other dim=5 
operators (in particular, with the clover term) is forbidden by the symmetries. 
The spectral quantities are therefore automatically 0{a) improved. 

The exact global symmetries above imply the existence of conserved cur- 
rents. The form of the conserved currents is not unique. It is very useful 
to work with conserved currents and scalar and pseudoscalar densities which 
transform covariantly (i.e. the same way as in the formal continuum) under the 
global transformations in Eqs. (|l9|) and (po[). These dim=3 operators are again 
automatically 0{a) improved: there are no dim=4 operators they can mix with. 

Our way to find chiral covariant conserved currents is the same as that 
proposed by Kikukawa and Yamada ]58| who presented explicitly the vector 
and axial currents in the overlap construction with Wilson kernel. We present 
here the currents in the general case in a form which, we believe, is easy to use 
in numerical simulations. We discuss the scalar and pseudoscalar densities, the 
Ward identities and the case 2R ^ 1 also. In the Appendices we collect some of 
the identities implied by the GW relation. 



5.1 A useful form of the currents 

Consider a global transformation ip — > ip + Sip, ip — > ip + Stp and assume that the 
action Eq. (|23|) is invariant under this transformation, 5 A = 0. Assume that 
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under the corresponding local transformation the change of the action can be 
written in the form 

SA = i^A(U)[D(U),e}B(U)ip , (24) 

where, for notational convenience, we treat the infinitesimal, x-dependent pa- 
rameter of the transformation as a diagonal matrix 

= ■ (25) 

We shall consider the case where A(U) and B(U) transform correctly under a 
gauge transformation to make 5 A gauge invariant. The corresponding current 
is defined through 

5A = iY,9A*)M*) = -*^<*KM*)> (26) 

X X 

where and d* are the forward and backward lattice derivatives respectively. 

Consider first the case A = B = 1, the generalization is trivial. We extend 
the gauge fields from SU(7V C ) to SU(iV c ) x U(l) maintaining the gauge covariance 
of D(U). Consider a U(l) gauge transformation 

U^x) - Ufj,(x) = e fc Wc^(a;)e- fc ( a+ « . (27) 

For this we have 

D(U) yz -> D(U) yz = e^D(U) yz e-^ , (28) 
i.e. for an infinitesimal U(l) gauge transformation the change of D(U) is given 

by 

5 s D(U) = -i[D(U),e\. (29) 
Using Eqs. @ and @ one has 

SA = i^d^{x)J^(x) = -^6 g D(U)ip. (30) 

X 

We show below that S g D(U) can be rewritten in the form 

8 g D(U) = -iY,9A*) K n(*)> ( 31 ) 

X 

which gives then 

J^x) = $Kp(x)1> = Yjv{K^x)) y ^ z . (32) 
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To construct K^x) satisfying Eq. (|3l]) observe that the gauge transformation 
in Eq. rt28) can be reached by performing the changes 



U^x) - U^(x) = e ia ^U^x) (33) 
on each link independently and taking the actual values of a^{x) to be 

atp(x) = -d^e{x) . (34) 



(Note that the individual changes in Eq. (]33j) are not pure gauge transformations 
- they only add up to that after all the links have been properly changed.) 



To linear order we have 



S e D(U) = D(U) - D(U) = -J2 M*) 



5D{U [ « ] ) 



a=0 



Eqs. (|3l|) and (|3j) give the kernel of the current 

K^x) = 



5D(ul a) ) 



Sa^(x) 



(35) 



(36) 



a=0 



It is easy to see that in the general case of Eq. (|2j 
the kernel 

Kp{x) = A(U)K^x;U)B(U) 



the current is given by 



(37) 



where K^{x) is defined in Eq. (|36|). Observe that Eq. (|36|) provides a straight- 
forward way for practical determination of the kernel K^(x) (and hence of the 
conserved currents discussed below) by performing the numerical differentiation 
in a M (x). 



5.2 Chiral covariant conserved vector and axial currents 

Consider a global chiral transformation acting only on the right-handed compo- 
nents 4>r = pRip and ip R — 4>Pl ■ 

5 a R ^ = ieT a P R ^, 5%~$ = -i$P L T a e (global), (38) 
and the analogous global left-handed transformations 

5l^ = ieT a P L ^j, 5 a L ^^-i^P R T a e, (global). (39) 
The projectors above are defined as M 



P R/L = -(1±75) 



P R/L = -(1± 75 ). 



(40) 
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These transformations are symmetries of the action. It is convenient to promote 
these global transformations to local ones in a way preserving chirality pifl , 

S a R ?P = iT a P R eP R ^ , 5<0 = ~i$P L eP L T a . (41) 

In Eq. (^d]) epRip is not a right-handed field if e is x-dependent. This explains the 
presence of the second projector P R Keeping the form of Eqs. (|38|), ( |39| ) in 
the local case would produce conserved currents as well, transforming, however, 
non-covariantly under global transformations. Note that it is preferable to use 
the covariant conserved current. For example, the pion is created from the 
vacuum by the covariant current, and the corresponding relation defining /„. 
is valid for this current. Non-covariant conserved currents result in an extra 
m-dependent factor (cf. Appendix C). 

Using the identity^] DPr = P^D one can write the corresponding change of 
the action as 

5%A = # (P L [D, e]P R T a ^j if, . (42) 

This is of the general form considered in Eq. (|24|) hence we readily obtain the 
corresponding current 

J a Ril (x) = ^K Rfl (x)1> , K% (x) = P L K, (x)P R T a , (43) 

where the kernel K^{x) is defined in Eq. (|36"|). 

Similarly, the left-handed local transformation is 

6 a L Tjj = iT a P L eP L ^ , 5$$ = -i$P R eP R T a . (44) 

and the left-handed current is given by 

Jl^x) = , K^{x) = P R K^x)P L T a . (45) 

Under the global chiral transformations in Eqs. ( p8| ) and (|3^) these currents 
transform covariantly 

~<^4» = ifab c J c Rfi (x) , -F L J h Lt Sx) = if abc J c L ,(x) . (46) 

The currents J Ril and Jl m are, of course, invariant under the left- and right- 
handed transformations, respectively. These properties imply that the vector 
and axial currents 

V;{x) = J£» + Jl^x) = i> (PlK^x)P r + P R K^x)P L ) T a V 

1_ ( 47 ) 
= -^(K^(x)- l5 K^x)%)T a i,, 

4 This and other useful identities are collected in the Appendices. 
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W = J hM - J W = 1> (PlK,(x)Pr - P R K,{x)P L ) T a tp 

1- 
2 1 



-V (-75^(^+^(^)75) TV, 



(48) 



are transformed under global vector and axial rotations 

S v =S R + 5 a L , S A = 6 R -S a L , (49) 

covariantly 

% -8 v V b Jx) = if a Mx) , -S a v A»(x) = if abc A c Jx) , 



e 



U a A V*{x) = ifabcA^x) , U a A Al{ X ) = ifabcV^x) 



(50) 



Let us discuss the overhead related to the currents in Eqs. ( ]47| ) and ( [48| ) 
in comparison with using the non-conserved, non-covariant local currents. The 
action of the kernel K^(x) on a vector v is oc D(U<- a >)v - D(U) v in the simplest 
numerical approximation of the derivative in Eq. (p6|). The projectors P R , Pl 
require an additional multiplication with D(U). We believe, this small overhead 
is justified by the advantages of having conserved and O(a) improved currents. 

Note that the fermion fields transform under the local vector and axial trans- 
formations as 

^ = ^ a i(e + 75 e 75 )V, 5^P = -i^eT\ (51) 
6 A iP = iT a ^{j 5 e + e%)iP, Sffi = i^e l5 T a . (52) 



5.3 Chiral covariant scalar and pseudoscalar densities 



It is easy to show that the scalar and pseudoscalar quantities 

S a = ^(l-^D\T a ^ : P a = v; 75 - ^) T<V , (53) 

a = 0,1,... Nj — 1 (with T° = 1) transform under global vector and axial 
rotations like in the formal continuum, 

l -8 a v S b = if abc S c , l -5 a v P b = if abc P c , (54) 



U a A S b = -d abc P c , U a A P b = -d abc S c , (55) 
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where {T a ,T b } = d a bcT c . In particular, for a flavour singlet axial transforma- 
tion we have 



-5 A S a 



-2P a . 



-5 A P a = -2S a : 



(56) 



while the non-singlet axial transformation of the flavour singlet quantities reads 



-5%S = -2P a , -6 A P=-2S a , a = l,...Nj-l. (57) 
(Here 6 a = 6 A , S — S°, etc. for the flavour singlet quantities.) 



Since the quantity S enters the action in the mass term mS, we need the vari- 
ation of S under a local chiral transformation when considering Ward identities. 
They read 



'5e{x) 



-2P a {x) 



. 8%P 
l Se(x) 



-2S a {x) , 



(58) 



where S a (x) and P a {x) are the covariant scalar and pseudoscalar densities 
related to the divergencies of conserved covariant currents (ff7|), (|4§|). Using 
Eqs. ©,@,(gg|) one obtains 



S a {x) = ip 



\e{ x ) (i - \d 



(59) 



P a (x) = ip 



\E(x) l5 (l - \D 



T a i\>. 



(60) 



Here we have introduced the notation 

( E ( X ))yz = 



(61) 



When summed over the lattice these densities reproduce the quantities in 
Eq. ©: 



(62) 



as it should be since 8/5e(x) corresponds to an infinitesimal global transfor- 
mation. 
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Although Eqs. (|59|) , (|60|) look inconveniently complicated, up to contact terms 
and prefactors (1 + (am)) they can be replaced by the corresponding point-like 
densities (cf. Appendix C). 

Note that under a global transformation the first and second terms in Eqs. (|5 
^0|) transform separately, i.e. the simpler expressions 

i>E{x) (l - ±Dj T a ^ , ^E(x) l5 (l - ±D) T> (63) 

also define covariant densities. Similarly, the non-conserved currents 

^E(x)^ (l - ±D) T a ^ , $E(x) W (l - ±D) I> (64) 

are also transforming covariantly. They are not related, however, to each other 
by a Ward identity. 



5.4 Ward identities 



Consider the fermion action with a flavour invariant mass termP] 

= A + mS = ^(D + m(l-~D S J S J-il). 

Under a local axial transformation defined in Eqs. and (49) we get 

i|^=9;4;(x)-2mP». (66) 

Consider the local Ward identity obtained by the change of variables defined by 
an axial transformation in the path integral for the fermionic expectation value 
of a multi- local operator 0(yi, 1/2, • • • , Vn)' 

( l W)) F ~ (°^ A >^F + 2m (° pa ( X ))F - *°°2JV) (Oq(x)) F = , (67) 

where P(x) is given by Eq. ( |60| ) and the un-normalized fermionic expectation 
value (in a given background gauge field) is defined as 

(0) F = [ [##]e-^0(^,V) ■ (68) 



3 In our convention the Boltzmann factor is exp(— Am)- 
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The last term in Eq. (|67j) is the contribution from the measure which is not 
invariant under a flavour singlet axial transformation and its change is 
given by 

[#'] -»• [d^ 1 ] = [#](!+ fy) (69) 

wher^J 

S(* = -iTi (PrcPr - P L (Pl) = -iTr (%e) = 2iN f e(x)q(x) , (70) 



or 



Se(x) 



-2N f q(x). (71) 



(Note that the [dip] part of the measure is invariant.) In Eq. ( |70| ) q{x) is the 
topological charge density 

q(x) = ^tr( 75 D{x,x)), (72) 
which enters the index theorem ||Tof 

tr(7 5 -D(a;, x)) = N f indcx(L>) . (73) 

X 

If x in Eq. ( |67| ) is sufficiently far removed from the operator 0(yi, y2, ■ ■ ■ , y n ) 
{\yi — x\ is much larger than the range of D), the first term in Eq. ( |67| ) is zero. In 
this case the Ward identity is consistent with the classical equations of motion 

d;A«(x) = 2mP a {x) (74) 

for a ^ 0. Summing over x in Eq. ( |67| ) leads to the global Ward identity 

.5°kO 



+ 2m{OP a ) F -6 a0 2vN f {O) F = 0, (75) 

F 



Sc 

where v is the value of the topological charge of the given gauge configuration. 



Let us illustrate the consequences of Eq. ( |75[ ) on two examples. Consider 
first a = 0, O = 1: 

rn{P) ¥ = vN f {l) F . (76) 
Combining this relation with that obtained by setting O = P: 

m(PP) F - (S) F - vN f {P) ¥ = . (77) 



6 To get the correct sign keep in mind that ip, ip are Grassmann variables. 
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we get 

v 2 Nj{l) F = m 2 {PP) F -m{S) F . (78) 

Averaging over the gauge fields leads to an identity for the topological suscep- 
tibility 



m 



m 2 1 



Xtop = Y = -N] {Sm + N]V {PP) ■ (?9) 

Assuming that there exists no massless excitation in the flavour singlet chancl 
we get in the chiral limit |30| 

(80, 

In the second example consider the non-singlet Ward identity with O = P b 

2m(P a P b )=d abc (S c ) , (81) 

and in particular 

2m(P 1 P 1 ) = j ff {S). (82) 

(Note that duo — V^/)- I n the chiral limit the pseudoscalar correlator is 
saturated by the Goldstone boson pole, while 

(S) = V(S{0)) = -VEN f (83) 

leading to the GMOR relation on the lattice |TTJ 

f 2 ml = 4mS . (84) 

5.5 The case of the general GW relation 2R 1 

Consider the Dirac operator satisfying the general GW relation 

75-D + D 75 = Dj 5 2RD . (85) 

In order to connect this with the 2R — ► 1 case we rescale D and the field 
variables: 

D = (2i?r 1/2 Di (2i?)~ 1/2 , (86) 

and 

^^ 1 {2R) 1/2 , = (2R) 1/2 ^ . (87) 
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Obviously, we have 

A(m) = %l> 1 

where 



Di+roj l--£>i 



ipi = ipD(m)ip . 



D(„,) = D + m \^-\d 



(88) 



(89) 



The new Dirac operator D\ satisfies the GW relation with 2R = 1 

75-Dx + Ci75 = #175^1 ■ (90) 

The spectrum of D\ lies on a circle of radius 1 with the centre at 1. Assuming 
that D is normalized in the standard way, so that for the free case its FT is 
D(p) = i-f^p^ + 0(ap 2 ), the operator D\ will be normalized differently: Di(p) = 
cijfiPn + C(ap 2 ) with c = 2R(0). This could be restored by a simple rescaling 
of the fields, however, for convenience we shall rather keep the non-conventional 
normalization of D\. Of course, this choice does not affect the physical results. 

One can derive currents and densities by writing them first in terms of ipi, 
tpi and D\. We choose here, however, a more direct way, and work with the 
original variables tj), ip and D. To start with, we introduce the operators 

f 5 = (2i?) 1 / 2 75 (2i?)- 1 / 2 = 75 (1 - 2RD) , (91) 

£fl = i(l + f s ), P £ = ~(l-f B ), 

which are related to their counterparts at 2R = 1 by a similarity transforma- 
tion. Although, these operators are not hermitian, their basic properties are 
unchanged: 

l5 D = -Df 5 , T 2 5 = l, VrVr^Vr, VrP l = 0, V R + V L = l. (92) 

It is convenient to define the local transformations through (cf. (|5l]), J5^)) 
5^ = iT a ^(e + het^, 8$$ = -i^eT a , (93) 

^ = l T a l (f 5 e + ef 5 ) V, = i^e l5 T a . (94) 

The corresponding global transformation is a symmetry for the massless case. 
The conserved currents are given by expressions analogous to Eqs. (ff7|), ( [f8| ) 
with the replacements Pr/l ^ Pr/Li 75 ~~ * 1*5, and where K^x) is also given by 
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Eq. (pq). It is straightforward to verify that they are transforming covariantly 
under global transformations. 



Similarly, the covariant densities obtained from 



are given by 

S a (x)=Jj 



T a i/j , 



(95) 



(96) 



P a (x) = i> 



T a i). 



(97) 



The presence of 2R ^ 1 leads to a small overhead only: 2i? and (2R)~ 1 are 
local, non-singular operators without Dirac indices. Since the inverse of D(m) 
can be written as 



Dim)- 1 



1 



1 - to/2 



2i? 



L>2i? 



1 - m/2 



(98) 



the multi-mass solver can be easily generalized to this case. The eigenvalue 
equation for Di leads to the generalized eigenvalue equation for D 



where the eigenvector \cj)\) is related to that of D\ by 

|0 A ) = (27?) 1 / 2 |^) , 
and forms an ortho- normalized system with the weight (2i?) _1 

{4>'\\^\<t>\) = S\\> ■ 



(99) 
(100) 
(101) 



Note finally that the presence of the extra factors appearing in eqs. (p6]),(p7|) 
in fact simplifies the calculation of expectation values because due to the identity 
(cf. 



— - -d] Dhn)- 1 = ] —- ; 

y 2R 2 J y ' (l-m/2) 2 



-1 -1 



D2R 



1 - m/2 



rn 



(102) 
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it removes the factor 2R in ( p8j ) . 

One can show that the covariant densities appearing in the Ward identities 
can be replaced also in this case (up to contact terms and am dependent factors 
going to 1 in the continuum limit) by simpler operators (cf. Appendix C) 

6 Conclusions 

The main purpose of this study has been to show that it is feasible to use the 
parametrized fixed-point QCD action in simulations. We also give a practical 
and general construction of conserved currents and covariant densities for chiral 
lattice actions. The tests of the parametrized fixed-point Dirac operator D FP 
show that the deviations from chiral symmetry are small and can be removed 
by small corrections in a straightforward fashion with the overlap construction. 
A first study of the hadron spectroscopy shows that the additive quark mass 
renormalization and, more importantly, the fluctuation of it, are small, allow- 
ing us to go quite small physical quark masses. The speed of light, extracted 
from the momentum-dependence of the light hadron spectrum, is consistent 
with 1, evidence that the fixed-point properties are intact. We also measure 
the renormalized chiral condensate directly using finite-volume scaling, giving 
£jyfg(2 GeV) = (262 ± 12 MeV) 3 , which is in good agreement with other recent 
measurements. In addition, we measure the quenched topological susceptibil- 
ity Xt = (196 ± 6 MeV) 4 and test other methods of determining S, using the 
pion mass measurements or measuring the distribution of the smallest non-zero 
eigenvalue of the Dirac operator. We also examine near-zero modes of the Dirac 
operator and find that they do appear to be chiral locally, as other studies have 
found, in support of the picture of instanton-dominancc. 

Any chiral lattice action is much more expensive to use in simulations than 
the standard actions. How competitive is the fixed-point QCD action with, say, 
the overlap operator? The majority of the simulation time is spent inverting the 
Dirac operator, so using the more expensive fixed-point gauge action 5*J P , with 
its desirable properties, is a small part of the overall cost. If one is interested 
in observables where small deviations from chiral symmetry are acceptable, the 
parametrized operator D FP can be used. An optimal implementation of D ov 
with D w in the kernel, say the rational approximation method, which is on 
the order of 100-200 times more costly than using Z? w , has chiral symmetry 
violations which are orders of magnitude smaller. The overlap- improved -D F ^T, 
constructed with a low order polynomial, achieves the same accuracy of chiral 
symmetry at a similar cost. Comparison of overlap and domain wall fermions 
does not show a large difference in cost for a given accuracy of chiral symmetry 
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pl| . The major advantage of the fixed-point action is if the lattice spacing 
dependence is much reduced. If this is so, the cost of using the fixed-point 
action is offset by being able to work on much coarser lattices. A large scale 
and systematic study, which is in progress, is required to accurately determine 
how large the lattice artifacts are. 

Quenched QCD simulations with chiral fermions have advanced very quickly 
over the last few years. However, the serious problem of how to implement chiral 
fermions in full QCD will have to be tackled. It is also an open question, how 
necessary chiral lattice fermions are for much of QCD phenomenology. 
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Appendix A. 



Some of the useful identities, which follow from the Ginsparg- Wilson relation 
are collected here (case 2R = 1). Concerning the definitions of Pr, Pi,, 75, Pr 
and P L we refer to Eqs. (||), Q. 



75D = -D% . 



(A.l) 



DP R = P l D, DP L = PrD 



L = ^Rl 



(A.2) 



D = P R DP L + P l DPr . 



(A.3) 



2P R P L = PrDP l = P R D = DP L , 
2PlPr = PlDPr = P L D = DPr , 



(A.4) 



PrPr = Pr\1-^D 



^-^d)p r . 



(A.5) 



PlPl = Pl I 1 - -D I = 



l-- 2 D)P L . 



(A.6) 



PrPr + PlPl =l--D, 



(A.7) 



PrPr - PlPl =^\\--D 



1 - 2^ H 5 • 



(A., 



34 



D{ »' ] = D " ( 1 - \ D ) • 1 \ D 1^72 (' - ■ (A - 9) 



Appendix B. 



We consider here identities which follow from the general GW relation with 
2R ^ 1 . Concerning the definitions of f^, Vr and Vl we refer to Eq. (BH). 



j 5 D = -DT 5 , (B.l) 



DV R = P l D, DT l = PrD , (B.2) 



D = P R LrP L + P L DP R , (B.3) 



^PrTl = PrDVl = P R D = DV L , 
Art 

^PlVr = P l DVr = P L D = DPr , 
ZH 









(±- 




(—--d) 


\2R 




\2R 2 J 



2R 













(—--d) 


\2R 




\2R 2 



(B.4) 



(B.5) 



1 PrPr =Pr[^-^D) = [^--^D)Vr, (B.6) 
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m p ^ = p -{m-\ D ) = {m-\ D ) p ^ (R7) 

^{PnV R + P L V L )=^- 1 -D, (B.8) 
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Appendix C. 



We show here, first for the case 2R = 1, that up to contact terms[] the covariant 
densities (p^), ( p0| ) can be replaced by their point-like counterparts. 



According to Eqs. (|60|), (|A.9|) one has 



P a {x) 



1 * 



1 - to/2 



-E{x) l5 1 - -D(m) 



1 



+ - 1--Z>(m) (%E(x) + E(x)%) 



(C.l) 



In the matrix elements one has to calculate D(m) 1 [. . . ]D(m) 1 and the terms 
proportional to D(m) can be omitted: 



1 - to/2 



(C.2) 



where the ~ sign means "equivalent up to contact terms" . The following rela- 
tions can be easily verified 



1 



75 



l-m/2 

1 + m/2 75 + l + m/2 



, . l + m/2 
D(to)7 5 = -75 



1 



l-m/2' l-m/2 



j 5 D(m). (C.3) 



For the %E(x) term in ( |C.2| ) we use the first relation while for E(x)j5 the second 
one. In this case one of the propagators is eliminated by the corresponding 
D(m), leading again to contact terms. The final result is 

pa (*) * ~ ^j^Eix^T^ 

(C.4) 



(1-to/2)(1-to 2 /4) 



1 



1 



ijE{x) lb 1 - -D T a ip . 



1 - TO 2 /4 

Note that these prefactors are 1 + O(etm) (restoring the lattice spacing a), i.e. 
even omitting them does not affect the continuum limit. 



The analogous calculation for the scalar density gives 

S a (x) « 1 _ 1 m/2 lpE(x)T a ip^E(x) (l-^D) T'r. (C.r„ 



Measuring the correlator of the covariant pseudoscalar density one can di- 
rectly determine the pion decay constant f n from the relation 

fl 3 

y<0|P a (f,t)P 6 (0,0)|0) = -h^-6 ab e~ m ^ . (C.6) 

16m z 

X 

1 To be understood in a broader sense: not necessarily proportional to S xy but negligible 
when \x — y\ is larger than the range of the Dirac operator. 



37 



(The normalization corresponds to the GMOR relation (p^), f%m„ = 4mS, i.e. 
U~130MeV.) 



Finally, we give the generalization of Eqs. ( p.4| ), (C.5) for the case 2R ^ 1: 
1 1- 



P a {x) 



l-m/2 2 



1 l + m 2 /4 1 
E{X) 2R + l-m 2 /4 2i? 



B(a:) 



75 T a V, (C.7) 



and 



1 1- 



l-m/2 2 



Sfa;)— + — -B(x) 
y '2R 2R { ' 



T a i\). 



(C.8) 
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1 2 

Figure 1: The full eigenvalue spectrum of D FP for a S*J P gauge configuration 
of volume 4 4 (crosses) and the smallest 100 eigenvalues on a 8 4 volume (circles) 
at a lattice spacing 0.16 fm. 
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Figure 2: Ratio of the 50 smallest A^A eigenvalues to the largest eigenvalue 
using D FF (filled squares) and £) Wllson (open squares) for 5 different S FP gauge 
configurations of volume 12 4 at a lattice spacing 0.16 fm. 
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Figure 3: Breaking of the Ginsparg- Wilson relation as measured by Aqw(-W) 
on 10 4 lattices at different lattice spacings and different number of exactly pro- 
jected A^A eigenvalues. Keeping the number of exactly projected A^A eigenval- 
ues constant, the exponential fall-off is steeper for smaller lattice spacings. If the 
number of exactly projected A^A eigenvalues is adjusted such that the range of 
the not exactly treated A^A eigenvalues is approximately the same for the dif- 
ferent lattice spacings, the fall-off is roughly equal in all cases, which illustrates 
that the approximation is governed by the ratio X m i n (A ji A) / X max (A ji A) . 
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Figure 4: The dependence of the deviation of the eigenvalues from the 
Ginsparg- Wilson circle on the order of the polynomial approximation, as mea- 
sured by |Re(A)|. 
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Figure 5: The locality of D ov , as measured by the expectation value f(r)/f(0) 
defined in Eq. (|), using D FP and D Wilson as input on 12 4 lattices at a = 0.16 fm. 
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Figure 6: ir, p and N masses for the parametrized FP Dirac operator D FP 
at (3 = 3.0 from 70 gauge configurations with lattice size 9 3 x 24. The dashed 
line shows a quadratic fit to the rho meson, leading to a lattice spacing of 
o~0.17£m. 
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Figure 7: Edinburgh plot for the parametrized FP Dirac operator on 9 3 x 24 
at [3 = 3.0 for the naive and covariant mass definition. 
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Figure 8: Squared pion mass for the parametrized FP Dirac operator D FP from 
100 6 3 x 16 gauge configurations. The inset shows the smallest four masses on 
a magnified scale to identify topological finite- volume effects. 
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Figure 9: Squared pion mass for the parametrized FP Dirac operator D FP 
from 70 9 3 x 24 gauge configurations. The massive Dirac operator is constructed 
with the covariant scalar density as in Eq. (88). Shown are quadratic fits to the 
P correlator at large and the P-S correlator at small quark mass with (dashed 
line) and without (solid line) the QxPT logarithm. The inset shows the smallest 
four masses on a magnified scale. Topological finite- volume effects are reduced 
considerably. Also plotted is the unrenormalized quark mass from the axial 
Ward identity (AWI) together with a linear fit. When the quenched chiral 
logarithm is included in the fit, the fits from the AWI quark mass and from the 
squared pion mass agree in the chiral limit. 
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Figure 10: Squared speed of light c 2 for pion and rho mesons with the 
parametrized FP Dirac operator D FP + m at the smallest non-zero momen- 
tum. The error bars are from statistical errors only. Systematic uncertainties 
from choosing the fit range (especially for the rho) would increase the errors 
considerably. 
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Figure 1 1 : Squared pion mass for the overlap- improved FP Dirac operator 
from 28 9 3 x 24 gauge configurations at lattice spacing 0.16 fm, with a quadratic 
fit to the P correlator at large and the P-S correlator at small quark mass. The 
inset shows the smallest four masses on a magnified scale. 
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Figure 12: Squared pion mass for the overlap- improved FP Dirac operator D, 
from 32 9 3 x 24 gauge configurations at lattice spacing 0.13 fm, with a quadratic 
fit to the P correlator at large and the P-S correlator at small quark mass. The 
inset shows the smallest four masses on a magnified scale. 
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Figure 13: The ratio (— (^^)m h vQ a3 )/( ma ) versus ma for different volumes 
and topological sectors. 
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Figure 15: Comparison of different measurements of rgSj^g(2 GeV) 



52 



10.0 



8.0 



6.0 



Q_ 



|Q|=0 fit 

- |Q|=1 fit 
|Q|=2 fit 




0.5 



Figure 16: Distributions of the smallest eigenvalue A of for topological 
sectors |Q| = 0,1,2, volumes 4 4 at lattice spacing 0.30 fm. The curves are the 
random matrix theory fit to the histograms. 




Figure 17: Distributions as in Fig. [l6| for volumes 4 4 at lattice spacing 0.22 fm. 



53 



o 

?5 



0.10 



0.08 



0.06 



0.04 



0.02 



0.00 



0.0 




■ FP 

ACHIR 

• CP-PACS (Iwasaki) 
o CP-PACS (LW) 

♦ SCRI 

t Boulder 



0.1 



0.2 



(a/r o r 



Figure 18: The quenched topological susceptibility measured by different ac- 
tions and techniques. The CP-PACS data are obtained by cooling using the 
Iwasaki and the Luscher-Weisz actions (see The other data are obtained 

by chirally symmetric actions. FP denotes the present determination using . 
The others are obtained using a chirally improved action Q(CHIR), the over- 
lap with the Wilson action |63)(SCRI) and the overlap with the planar action 
@ H (Boulder). 



■54 



0.10 



0.08 



10% 

5% 

1% 



0.06 



x 

Q_ 



0.04 



r 1 -r ' 



0.02 




0.00 

-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

X 

Figure 19: The distribution P{X) for volumes 10 4 at lattice spacing 0.13 fm. 
We include the top 1%, 5% and 10% lattice sites with largest density ^^(x). 
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